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A CHANDRA OBSERVATION OF M51: ACTIVE NUCLEUS AND NUCLEAR OUTFLOWS 
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ABSTRACT 

We present a Chandra ACIS-S observation of the nuclear region of the nearby spiral galaxy M51 (NGC 
5194), which has a low-luminosity Scyfcrt 2 nucleus. The X-ray image shows the nucleus, southern 
extranuclear cloud, and northern loop, the morphology of the extended emission being very similar to 
those seen in radio continuum and optical emission line images. The X-ray spectrum of the nucleus is 
well represented by a model consisting of soft thermal plasma with kT « 0.5 keV, a very hard continuum, 
and an Fe Ka emission line at 6.45 keV with an equivalent width of > 2 keV. The very strong Fe line 
and the flat continuum indicate that the nucleus is obscured by a column density in excess of 10^** cm~^ 
and the spectrum is dominated by reflected emission from cold matter near the nucleus. The X-ray 
spectra of the extranuclear clouds are well fitted by a thermal plasma model with kT « 0.5 keV. This 
spectral shape and morphology strongly suggest that the clouds are shock heated by the bi-polar outflow 
from the nucleus. The shock velocities of the extranuclear cloud and northern loop inferred from the 
temperatures of the X-ray gas are 690 km s"^ and 660 km s"^, respectively. By assuming a steady-state 
situation in which the emission of the extranuclear clouds is powered by the jets, the mechanical energy 
in the jets is found to be comparable to the bolometric luminosity of the nucleus. 



Subject headings: galaxies: active — galaxies: individual (M51) 
— X-rays: galaxies — shock waves 



galaxies: nuclei — galaxies: Seyfert 



1. INTRODUCTION 

M51 (NGC 5194) is a nearby (distance 8.4 Mpc; Feld- 
meier et al. 1997) spiral galaxy hosting a low-luminosity 
active galactic nucleus (AGN). The nucleus shows optical 
emission lines and is classified as a Seyfert 2 (e.g., Ford 
et al. 1985; Ho, Filippenko, & Sargent 1997). Radio con- 
tinuum and optical emission-line observations reveal a bi- 
polar structure comprising a southern extranuclear cloud 
(hereafter XNC) and a northern loop; this morphology and 
the high gas velocities observed are indicative of nuclear 
outflow (Ford et al. 1985, Cecil 1988, Crane & van der 
Hulst 1992). A radio jet connecting the nucleus and the 
XNC is also seen in a 6 cm image (Crane & van der Hulst 
1992). The optical emission-line clouds are inferred to be 
heated by the jet emanating from the nucleus. 

X-ray observations are a powerful tool to investigate the 
physical nature of these phenomena. However, the spatial 
resolution of previous observations is not good enough to 
resolve the structures in the nuclear region. In the soft X- 
ray regime, extended nuclear emission was detected with 
the Einstein HRI (Palumbo et al. 1985), the ROSAT 
PSPC (Marston et al. 1995), and the ROSAT RRl (Ehle, 
Pietsch, & Beck 1995; Roberts & Warwick 2000), but only 
a hint of enhanced emission in the direction towards the 
XNC is seen in the ROSAT RRl data. The spectrum of 
the nuclear region obtained with the ROSATPSPC is very 
soft. Hard X-ray emission, with an X-ray luminosity of 
9x10'*'^ ergs s~^ (2-20 keV) was detected by Ginga and 
interpreted as a combination of low-mass X-ray binaries in 
the host galaxy and a highly absorbed AGN (Makishima 
et al. 1990). Subsequently, ASCA detected a strong Fe Ka 



emission line with an equivalent width (EW) of w 1 keV, 
which probably comes from the heavily obscured AGN, 
as well as emission from thermal plasma with a tempera- 
ture of RiO.5 keV. The 2-10 keV luminosity of the entire 
galaxy measured by ASCA was 8.4 x 10'^^ ergs s~^ which 
is about six times smaller than that obtained with Ginga 
in the same band (Terashima et al. 1998). The reason for 
this large difference in luminosities has recently been clar- 
ified by BeppoSAX observations. Fukazawa et al. (2001) 
find that the nucleus is photoelectrically absorbed below 
10 keV, but seen directly through the absorption at 20 
keV, implying an absorbing column of Nh — 5.6 x lO^'' 
cm~^. The 2-10 keV luminosity measured by BeppoSAX 
was similar to that of ASCA. Fukazawa et al. (2001) at- 
tributed the higher Ginga luminosity to variability of the 
absorbing column, as long as the Ginga measurement is 
not contaminated by nearby bright sources. None of these 
X-ray studies has sufficient spatial resolution to resolve the 
nucleus and outflows seen in radio continuum and optical 
[N II]AA6548,6583 images. 

High spatial resolution X-ray observations with the 
Chandra X-ray observatory can easily resolve the multiple 
components in the nuclear region and enable us to study 
the physical nature of the nucleus and outflows, as well 
as the various X-ray sources in the host galaxy. In this 
Letter, we present Chandra results on the nuclear region 
of NGC 5194. Results on the diffuse X-ray emission and 
discrete X-ray sources seen in the host galaxy will be pre- 
sented in a future paper. At the distance of 8.4 Mpc, 1 
arcsec corresponds to 40.7 pc. 
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2. OBSERVATION AND DATA ANALYSIS 

M51 was observed on 2000 June 20 with the Chandra 
Advanced CCD Imaging Spectrometer (ACIS: Garmire et 
al. 2001, in preparation). M51 was located on the back il- 
luminated CCD chip S3. The data were reduced using the 
CI AO software version 1.1.5. Gain correction was made 
using the latest calibration files. We used only ASCA event 
grades 0, 2, 3, 4, and 6 (see Chandra Proposer's Observa- 
tory Guide, 2000). The background was stable during the 
observation and a net exposure time of 14865 sec was ob- 
tained. The count rates of all the detected sources in the 
field were relatively low and effects of pileup were negli- 
gible. Spectral fits were performed using XSPEC version 
11.0 and all the quoted errors are 90% confidence range 
for one parameter of interest. 

3. RESULTS AND DISCUSSION 

3.1. X-ray Images 

Our Chandra image of NGC 5194 shows X-rays from a 

bright nuclear region, diffuse emission distributed over the 
host galaxy, and many discrete sources. The companion 
galaxy NGC 5195 is also detected in the same field. In 
this Letter, we concentrate on the nuclear region of NGC 
5194. 

An expanded image of the central region of NGC 5194 
in the 0.5-8 keV band is shown in Figs 1 and 2. In Fig. 
1, red, green, and blue colors correspond to the energy 
bands 0.5-1.3, 1.3-3, and 3-8 keV, respectively. The nu- 
cleus, extended features, and a few discrete sources are 
seen. The extended structures north and south of the nu- 
cleus are remarkably similar to the features seen in radio 
(6 cm , 20 cm) and [N II] images, and we refer to them 
as the northern loop and the XNC, respectively. Fig. 2 
shows a 20 cm image (contours, from Ford et al. 1985) 
overlaid on a Chandra image in the 0.5 -8 keV band. Both 
images have « 1" resolution and a close correspondence is 
seen between the two maps. The nuclei were aligned by 
shifting the X-ray image by 0.4". 

The color (^hardness) map (Fig. 1) shows that the nu- 
cleus has a hard spectrum while the XNC and the northern 
loop are soft. A spectral analysis (next subsection) shows 
that the nucleus emits both hard and soft spectral com- 
ponents, the latter being very similar to the spectra of 
the XNC and northern loop. The north-east part of the 
XNC is somewhat harder than the rest of the XNC and 
the northern loop. It is intriguing to note that this re- 
gion coincides with a bright knot in the [N II] image from 
which the broadest optical emission lines (extending over 
~ 1500 km s^^) are found (the point P in Cecil 1988). 
It is also worth noting that the elongation of the nucleus 
seen in [N II] is towards the point P, and that the radio 
jet connects the nucleus and point P (Crane & van der 
Hulst 1992). Thus this medium-hardness region might be 
the site where the jet emanating from the nucleus is cur- 
rently interacting with the interstellar medium and shock 
heating it. 

3.2. X-ray spectra 

3.2.1. Nucleus 

The spectrum of the nucleus was extracted from a cir- 
cular region with a radius of 1.5 arcsec. The spectrum was 



binned so that each bin contains at least three counts in 
order to enhance the visibility of the narrow Fe-Ka emis- 
sion line. In this chi-squared fit is not adequate 
and a maximum-likelihood method using the C-statistic 
(Cash 1979) was employed. In the fit with the C-statistic, 
background cannot be subtracted, so we first constructed 
a model of the background spectrum. This background 
model was obtained from an area away from the nucleus, 
which includes diffuse emission in the host galaxy. The 
spectrum from this area was fitted by a two temperature 
(A:T=0.23 keV, 0.60 keV; abundance 0.16 solar) plasma 
model plus a flat power law. This fixed background model 
was then added to the model for the nucleus after normal- 
izing by the geometrical area. 

The spectrum of the nucleus (Fig. 3a) consists of at 
least three components, namely a soft component show- 
ing emission lines from ionized species, a hard component 
which dominates the 3 8 keV band, and a strong emission 
line at 6.4 keV. We modeled the soft component with the 
optically-thin thermal plasma model MEKAL. Two mod- 
els were examined for the hard component: an absorbed 
power law or pure reflection from almost face-on cold mat- 
ter (Magdziarz & Zdziarski 1995). The latter corresponds 
to a situation in which the X-ray emission from the nu- 
cleus is almost completely obscured by a column density in 
excess of 10^'* cm~^ and only emission scattered by mat- 
ter near the nucleus is observed. We added a Gaussian 
component to represent the emission line at 6.4 keV. The 
absorbed power law and pure reflection model provide a 
similar quality of fit (C-statistic = 39.1 and 39.6, respec- 
tively) . The spectral parameters obtained are summarized 
in Table 1. 

In the case; of the best-fit absorbed power law model, 
we fixed the photon index of the power law component 
at 2.0 (typical of unobscured AGN) and obtained an ab- 
sorption column density of A^H=3.lt}'5 x lO^'' cm^^. The 
emission line centroid energy and equivalent width (EW) 
are 6.458;^ !^!^^ keV (source rest frame) and 3.5^^;^ keV, 
respectively. Although this model can describe the data, 
the very large EW of the Fc line is not compatible with 
the relatively small absorption column (3 x lO^'' cm^^), 
from which an Fe line with much lower EW (^ 200 eV) 
is expected (e.g, Awaki et al. 1991; Leahy & Creighton 
1993; Ghiseinni, Haardt, & Matt 1994 ). 

The center energy, line width, and EW of the strong 
emission line in the pure reflection model are 6.458^0 048 
keV at the source rest frame, a = 0.020 (<0.099) keV, 
and 4.8^2 5 keV, respectively. This line centroid energy is 
consistent with fluorescent iron emission from neutral or 
low-ionization state (< FeXIX) iron when possible small 
calibration uncertainties are taken into account. The line 
intensity of i^A't'^'t) x 10~^ photons cm~^ is consistent 
with that obtained with an ASCA observation (Terashima 
et al. 1998), although the error ranges are large. The 
huge EW (>2.4 keV) is difficult to explain even in the 
pure refiection picture and probably an overabundance of 
iron is required (Matt, Fabian, & Reynolds 1995). Such 
a large equivalent width is also seen in other Compton 
thick Seyfert 2s such as the Circinus galaxy (Matt et al. 
1996). The X-ray luminosity of the reflection component 
in M51 is 9.8 x 10^^ ergs s"^ in the 2-10 keV band. Al- 
though this luminosity is nine times smaller than the total 
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2-10 kcV luminosity of M51 observed by ASCA, X-ray 
binaries contribute significantly to the ASCA luminosity 
(Terashima et al. 1998). Therefore, the smaller luminos- 
ity obtained with Chandra does not imply time variability 
of the AGN. Recently BeppoSAX observations have shown 
that the nucleus is heavily absorbed by a column density 
of A^H = 5.6 X 10^^ cm~^ and the intrinsic luminosity cor- 
rected for the absorption is 1.1 x 10^^ ergs s~^ in the 2-10 
keV band (Fukazawa et al. 2001). Thus, the luminosity 
obtained with Chandra is about 1% of the intrinsic; X-ray 
luminosity, implying that ~1% of the nuclear luminosity 
is scattered in our direction. 

The soft component (Fig. 3a) is well represented by a 
MEKAL plasma model (Table 2). The obtained absorp- 
tion column density is only slightly larger than the Galac- 
tic value (1.3 x 10^° cm'^. Stark et al. 1992) given the 
large errors. The emission measure is 1.2 x 10^ cm~^ and 
the luminosity corrected for the absorption is 4.5 x 10'^^ 
ergs s~^ in the 0.5-4 keV band. The temperature and 
abundance are quite similar to the extended emission (the 
XNC and north loop). The gas may be shock heated as 
discussed below. 

3.2.2. Southern Extranuclear Cloud and Northern Loop 

We extracted X-ray spectra of the southern extranuclear 
cloud (XNC) and the northern loop. A background spec- 
trum was obtained from a region free of the nuclear emis- 
sion. Bright point sources in this region were excluded. 
This background spectrum was then subtracted from the 
spectra of the XNC and northern loop. The spectra were 
binned so that at least 15 counts were contained in each 
bin, so chi-squared statistics can be used. 

The two spectra (Fig. 3b and c) show emission lines 
which are indicative of optically thin thermal plasma emis- 
sion with a temperature below 1 keV. Therefore, we fitted 
the spectra with the MEKAL model modified by photo- 
electric absorption along the line of sight. This model 
provides successful descriptions of both the XNC and the 
northern loop. The results of the spectral modeling are 
shown in Table 2. 

The spectral parameters obtained for the XNC are 
kT = 0.58to°ot keV, abundance 0.141°;°^ solar, and A^h 
= 9.5I3 9 X 10^° cm~^. We note that the abundance val- 
ues in all the fits in this paper are unusually low for the 
central regions of spiral galaxies and arc similar to other 
X-ray measurements of abundance of spiral galaxies and 
starburst galaxies. The reason for these unphysically low 
abundances is unknown, but may involve systematic errors 
in modeling thermal plasmas, nonthermal continuum con- 
tributions, or depletion of gas phase elements onto grains 
(Netzer & Turner 1997, Arimoto et al. 1997, Dahlem et 
al. 2000, Matsushita et al. 2000) and does not necessarily 
mean real low abundance. The luminosity (corrected for 
absorption) and emission measure are 1.5 x 10^^ ergs s~^ in 
the 0.5 4 keV band and 2.7 x 10^^ cm~^, respectively. If we 
assume a simple spherical shape with a radius of 80 pc and 
a filling factor r]{< 1) of the X-ray emitting gas, we can de- 
rive some physical quantities. The density n, gas mass M, 
and the thermal energy U are determined as n = 2.1r7~"'^/^ 
cm-3, M = 1.5 X 10^t]^/'^Mq, and U = 1.2x lO^^rj^/'^ ergs, 
respectively. 

This thermal plasma is expected from shock heating 



due to outflows. The expected postshock temperature is 
Tg = 1.4x 10^ K, where Vs7 is the shock velocity in units 
of 100 km s~^, if the preshock medium is fully ionized and 
the helium abundance is 10% of hydrogen by number (Hol- 
Icnbach and McKec 1979). Using this relation, the shock 
velocity is determined as 690 km s~^ from the observed 
temperature of the XNC. This shock velocity is roughly 
consistent with that determined from [N II] observations 
(500 km s-i; Cecil 1988). The XNC shows LINER-like 
optical emission lines (Ford et al. 1985). The emission 
line ratios [O I]/Ha, [O III]/H/3, and [S II]/Ha are in the 
range expected from fast shocks without precursors for a 
shock velocity of 300-500 km s"-*^, as calculated by Dopita 
& Sutherland (1995, 1996). The observed [N II]/Ha ratio 
is larger than predicted by the shock model and could be 
due to an overabundance of nitrogen, as often suggested 
previously. Thus the optical emission line ratios are consis- 
tent with fast shocks and reinforce a shock interpretation 
of the origin of the extranuclear clouds. 

The spectral parameters for the northern loop are 
kT = 0.53l[!;j]^ kcV, abundance 0.14tE5;55^ solar, and A^h 
= 6.2l3'3 x lO'^'^ cm~^. The luminosity (corrected for ab- 
sorption) and emission measure are 2.0 x 10^® ergs s~^ in 
the 0.5-4 keV band and 3.8 x 10^^ cm~^, respectively. As 
in the case of the XNC, we determine physical parameters 
of the X-ray gas by assuming a simple geometry. Here we 
assume a spherical shell with an outer radius of 200 pc 
and a thickness of 40 pc, and a volume filling factor of ri. 
The parameters obtained for the hot gas are n = 0.88rj~^''^ 
cm-3, M = 5.0 X lO^ryi/^Mg, and U ^ 3.7 x IO^^t]^/'^ ergs. 

The X-ray spectral shape (i.e. the model's kT and abun- 
dance). X-ray morphology, and optical emission line prop- 
erties of the northern loop are quite similar to the XNC 
and suggest a similar origin to the XNC, i.e., shock heating 
by nuclear outflows. The shock velocity required to heat 
the gas to kT =0.53 keV is 660 km s~^, almost the same 
shock velocity as obtained for the XNC. The fact that the 
northern loop is much larger than the XNC suggests that 
it is an older structure. 

The total radiative luminosity of the shock heated gas 
(I/t) can be estimated from the H/? luminosities (Ford et 
al. 1985) by assuming the shock model (Dopita & Suther- 
land 1995, 1996). A value Lt ~ 2 x 10^^ gj.gs g-i jg 
obtained. This luminosity is comparable to the photon 
luminosity of the nucleus if an intrinsic X-ray luminosity 
of 10*^ ergs s~^ (Fukazawa et al. 2001) and a bolomet- 
ric correction of a factor of ten (Ho 1999) are assumed. 
If the shock-heated gas is, indeed, powered by the jets, 
which are assumed to be steady, the mechanical and ra- 
diative outputs of the nucleus are then similar to within 
the considerable uncertainties. 

4. CONCLUDING REMARKS 

Our Chandra observations have shown that the X-ray 
morphology of the extranuclear regions of NGC 5194 is 
very similar to that seen in radio continuum and opti- 
cal emission-line imaging. The X-ray spectra of the ex- 
tranuclear cloud and the northern loop are well described 
by thermal plasma models with kT ps 0.5-0.6 keV. This 
temperature is close to that expected from shock heating 
given the range of velocities seen (~500 km s~^) in opti- 
cal emission-line spectra. The optical emission-line ratios 
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also agree with models involving shocks with velocities in 
the range 300-500 km s~^. Further, the location with the 
hardest X-ray spectrum in the extranuclear cloud coin- 
cides with the region exhibiting the broadest optical lines 
1500 km s^"'^). This location is also apparently the ter- 
minus of the radio jet. All of these findings indicate that 
the extranuclear gas is shock heated by nuclear jets and/or 
other coUimated outflows. 

The nucleus itself exhibits a soft X-ray spectrum which 
is very similar to that of the extranuclear gas. We pre- 
sume that this gas is also shock heated by mass outflows. 
In addition, the nucleus shows a very hard continuum and 
an iron line with equivalent width in excess of 2 keV. This 



spectrum arises through reflection of the radiation from a 
compact nuclear source with a power law spectrum by cold 
matter in its vicinity. The nuclear source itself is obscured 
by a large column density of gas (7Vh=5.6 x 10^^ cm~^) 
and is only seen directly above 10 keV (Fukazawa et al. 
2001). 

We are grateful to Gerald Cecil for providing the radio 
and [N II] images in computer readable format. YT was 
supported in part by a Japan Society for the Promotion 
of Science Postdoctoral Fellowship for Research Abroad. 
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Table 1 

Spectral modeling of the hard component of the nucleus. 



Model 


[10^2 cm-2] 


photon index 


rpa 

[kcV] 


Width 
[kcV] 


EW 
[kcV] 


Luminosity'' 
[10^9 ergs s-i] 


C-statistic 


1 
2 


oi+i8 

0.25(< 34) 


2.0^ 
2.0^= 


D.40S_Q 

fi 4c:o+0.035 

D.^OO_o.o48 


0(< 0.085) 
0.020(< 0.099) 




1.1 
0.98 


39.1 
39.5 



Note. — Model 1: MEKAL + absorbed power law + Gaussian line, Model 2; MEKAL + pure reflection + Gaussian line, a: Source rest 
frame. 6: Observed luminosity (not corrected for absorption) in the 2-10 keV band, c: Fixed parameter. 



Table 2 

Spectral modeling of the soft thermal component. 



iVn kT abundance Luminosity" x^/dof 

[10^2 cm-^] [kcV] [solar] [10^9 ergs s"^] 

Nucleus" O.M^o'i; 0.49;o ^8 0-084 ± 0.05 05 Z 

XNC 0.095t|]'J5^i^ 0.581°-°^ 0.141°;°^ 1.5 53.9/47 

Northern Loop 0.062l°:g27 0.53lg:g^ 0-14to:gI 2.0 44.9/50 



Note. — a: Luminosity corrected for absorption in the 0.5-4 keV band, b: The hard component is modeled by a pure reflection continuum 
+ Gaussian line. 
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13" 29™ 54^ 53^ 52^ 51^ 



Fig. 1. — Chandra ACIS-S image of the central region of NGC 5194. Red, green, and blue colors correspond to the 0.5—1.3, 1.3-3, and 
3-8 keV band, respectively. The color intensity is scaled logarithmically to the number of photons in each pixel. The position of the optical 
nucleus is {13''29™52=.7, 47°11'43") (J2000). The size of the image is 40 arcsec X 45 arcsec. North is up and east is to the left. 




13'^29"54" 53"^ 52^ Si'' 

Fig. 2. — A VLA 20 cm map (contours, from Ford et al. 1985) overlaid on a Chandra image in the 0.5—8 keV band. Contours are plotted 
at 7.5, 15, 21, 30, 43, 60, 85% of peak and the peak corresponds to 2.12 mjy beam~^. 
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Fig. 3. — X-ray spectra of the nucleus (o), southern extranuclear cloud (XNC) (6), and northern loop (c). Crosses are data and histograms 
are models. The hard component of the nucleus is assumed to be reflected emission from cold matter (dotted line), while the soft component 
of all the three spectra is assumed to be thermal emission from a gas in coUisional-ionization equilibrium, modeled by the MEKAL code. 



